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Abstract 

Following a previous calculation of quark scattering in eikonal approximation, this paper presents 
a new, analytic and rigorous approach to the calculation of QCD phenomena. In this formulation 
a basic distinction between the conventional "idealistic" description of QCD and a more "real- 
istic" description is brought into focus by a non-perturbative, gauge-invariant evaluation of the 
Schwinger solution for the QCD generating functional in terms of the exact Fradkin representa- 
tions of the Green's functional Gc{x,y\A) and the vacuum functional L[A]. Because quarks exist 
asymptotically only in bound states, their transverse coordinates can never be measured with ar- 
bitrary precision; the non-perturbative neglect of this statement leads to obstructions that are 
easily corrected by invoking in the basic Lagrangian a probability amplitude which describes such 
transverse imprecision. 

The second result of this non-perturbative analysis is the appearance of a new and simplifying 
output called "Effective Locality", in which the interactions between quarks by the exchange of 
a "gluon bundle" - which "bundle" contains an infinite number of gluons, including cubic and 
quartic gluon interactions - display an exact locality property that reduces the several functional 
integrals of the formulation down to a set of ordinary integrals. It should be emphasized that 
"non-perturbative" here refers to the effective summation of all gluons between a pair of quark 
lines - which may be the same quark line, as in a self-energy graph ~ but does not (yet) include a 
summation over all closed-quark loops which are tied by gluon-bundle exchange to the rest of the 
"Bundle Diagram" . As an example of the power of these methods we offer as a first analytic calcu- 
lation the quark-antiquark binding potential of a pion, and the corresponding three-quark binding 
potential of a nucleon, obtained in a simple way from relevant eikonal scattering approximations. 
A second calculation, analytic, non-perturbative and gauge- invariant, of a nucleon-nucleon binding 
potential to form a model deuteron, will appear separately. 
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I. INTRODUCTION 



' The strong interactions comprise a richer field than the set of phenomena that we have 
learned to describe in terms of perturbative QCD or the (near-) static non-perturbative do- 
main of lattice QCD ( ■■■) It may well be that interesting unusual occurrences happen outside 
the framework of perturbative QCD-happen in some collective, or intrinsically nonperturba- 
tive way.^ It is in these terms that the Resource Letter: Quantum Chromodynamics, arXiv: 
1002.5032v2 [hep-ph], of February 26, 2010, concludes its overall review of the most salient 
achievements realized in QCD. Of course by its very definition, the non-perturbative qualifier 
applies to a large variety of very different realizations: There are several ways of 'being non- 
perturbative', in the either context of Lagrangian Quantum Field Theory or of Algebraic 
Quantum Field Theory, the former admittedly more pragmatic. 

It is this route of Lagrangian Quantum Field Theory which will here be followed so as to 
introduce an extraordinary property of fermionic QCD amplitudes. This property can be 
phrased as follows: 

In any Quark/Quark (or Anti-Quark) scattering amplitude, the full gauge-invariant sum 
of all cubic and quartic vectorial gluonic interactions, fermionic loops included, results in a 
local contact-type interaction mediated by a tensor field which, in both internal (color) and 
external (Lorentz) indices is rank 2 and antisymmetric. 

This property, dubbed 'Effective Locality' appears to be a genuine one of QCD and is 
worth exploring in relation to the expected non-perturbative properties of QCD, and some 
examples will be presented. It is most interesting that a previous approach, looking for a 
'most dual' description of QCD, and first restricted to the pure Yang-Mills case, was able to 
display some positive results in this direction [1]. It may also be worth noting that a non- 
perturbative analysis of Yang-Mills thermodynamics has recently put forth the occurrence of 
a contact-type interaction within the confined phase of SU{2) and SU{3) theories {2]. And 
even more recently, a fascinating functional connection of our approach to a most general 
formulation of quantum field theories [sl has been brought to our attention. 

In a previous paper it was shown how the formal Schwinger solution to the 

QCD generating functional can, by a slight rearrangement, be brought into gauge-invariant 
form. This simple procedure has been overlooked for decades. When combined with the 
Fradkin representations for the functionals Gc[A] and L[A] (that are an integral part 
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of the Schwinger solution), as well as with the useful Halpern representation jo], [fol for 
exp [— I / F^], the relevant Gaussian functional operations may be performed exactly. This 
corresponds to the summation of all Feynman graphs of gluons exchanged between quarks, 
and one then explicitly sees the cancelation of all gauge- dependent gluon propagators. Gauge 
invariance is achieved by gauge independence, a hope some authors long had {e.g. R.P. 
Feynman) for QED ll|, which turns out to be realized in the case of QCD thanks to its 
very non-abelian structure. 

One also can see the appearance of a new and exact property of Effective Locality (EL), 
which simplifies all calculations by transforming the remaining functional integrals into sets 
of ordinary integrals, a non-trivial mathematical point jsj. 

Further, it readily becomes apparent that one cannot continue to consider quarks in the 
same, conventional fashion as quanta of other {e.g., abelian) fields, such as electrons, which 
satisfy the standard measurement properties of quantum mechanics, perfect position de- 
pendence at the cost of unknown momenta, and vice-versa. This is impossible for quarks 
since they always appear asymptotically in bound states, and their transverse coordinates 
can never, in principle, be exactly measured. In the developments to follow one finds that 
a violation of this principle produces absurdities in the exact evaluation of any QCD am- 
plitude. At a phenomenological level, we therefore propose a change in the basic QCD 
Lagrangian, which introduces a measure of transverse fluctuations by the insertion of an 
unknown probability amplitude that is essential to and determined from quark binding into 
hadrons. All previous absurdities in estimates of all "realistic" QCD amplitudes are then 
removed, and one finds a simple method of constructing quark-binding potentials, as well 
as other, somewhat more involved potentials between nucleons to form nuclei. 

This paper is divided into a sequence of Sections each emphasizing a particular aspect of 
this formulation. In Section |Tll as a result of the Gaussian nature of the functional operations 
which are performed exactly, the expression of the fundamental gauge invariance becomes 
evident, along with the EL property and its simplifying effects. In Sections IIIIl IIVI and |Vl 
the obvious need for transverse fluctuations is made clear, and the form of a corresponding 
insertion of an unknown probability amplitude into the basic QCD Lagrangian is described. 
Sections IVTI IVIIt and IVIIII deal with the construction of quark-binding potentials, exhibited 
for a model pion and model nucleon, using a simple, analytic technique (which reverses the 
Potential-Theory derivation of an eikonal function from a specified potential). Section HXl 
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is devoted to a brief Summary and Speculations. In the remainder of the present Section, 
we shall for completeness summarize the argument of ^ , which shows how the conventional 
Schwinger solution for the generating functional may be converted, for QCD, but not for 
QED, into a manifestly gauge-invariant form. 

It is simplest to begin with QED, and its free-photon Lagrangian, 

C-o = -^^^,u^^,u = {df,A^ - d^A^f , (1) 
whose Action Integral may be rewritten as 



d'xCo = -- I {d.A,Y + - I {d,A,Y (2) 



The difficulty of maintaining both manifest gauge invariance (MGI) and manifest Lorentz 
covariance (MLC) appears at this stage. What has typically been done since the original 
days of Fermi, who simply neglected the inconvenient {d^A^Y , is to use the latter to define 
a relativistic gauge in which all calculations maintain MLC, while relying upon strict charge 
conservation to maintain an effective gauge invariance of the theory. 

The choice of relativistic gauge can be arranged in various ways. The simplest functional 
way is to multiply the inconvenient term by the real parameter A, and treat it as an interac- 
tion term. For definiteness, one can begin with the free-field, (A = 0, Feynman) propagator 
Dc°^y = Dc, where (—9^) Dc = 1; then the free-field Generating Functional (GF) is given 
by 

3r{j} = exp{^|j-Df.jj. (3) 
Then, operating upon it by the 'interaction' A-term, a new free-field GF is produced 

3f){j} = ei^/^^'^^")' .e^r^'"'-^ (4) 



where Di^Ji^ = [S^u — C^^idy/d'^] D^, with ^ = A/(l — A). The functional operation of (jl]) is 
fully equivalent to a bosonic, gaussian functional integration (FI). Such 'linkage operation' 
statements are frequently more convenient than the standard FI representations, since they 
do not require specification of infinite normalization constants. 
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The Tr-Log term of (jl]) is an infinite phase factor, representing the sum of the vacuum 
energies generated by longitudinal and time-like photons, with a weight A arbitrarily inserted; 
this quantity can be removed by an appropriate version of normal ordering, but can more 
simply be absorbed into an overall normalization constant. 

Including the conventional fermion interaction, £int = —igip'y ■ Aip, and the gauge 'inter- 
action' (d^A^)'^, it is easy to show 0, Q] that the standard Schwinger solution for the GF 
in the ^-covariant gauge becomes 

3g)^[^-,,^,,^] = A/-e*/^-Gc[A]..+L[A]+iA/(a,A,)^ .el/^-D^"'-^, (5) 

A— 

where Gc[A] = [m + 7 ■ (9 - igA)]-\ L[A] = Tr In [1 - ig-f ■ ASc], Sc = Gc[0], and where 
the phase factor of (jlj) has been absorbed into A/". A most convenient re- arrangement of ([5]) 
uses the easily-proven identity [7], for an arbitrary functional ^[^], 



_i 5j 

where Ti'f' = -| / ^ ■ D^^^ ■ ^, so that now reads 



A=J T>i>'-j 

This is the formal solution for the GF of QED in the gauge ( that has been known and used 
for a half-century j^]. 

We now come to QCD, with 

>Cqcd = -\f%F;, -^.[m + ^^{d,- zgA^X^)] ■ ^, (8) 

where F;^^ = d^^A^ - d^A"^ + gf^^A^A^ = f^^ + gf^^A^Al. One can rely on the following 
observations: 'proper' quantization in the Coulomb gauge, for the free and interacting theo- 
ries, yields the same equal-time commutation relations (ETCRs) for QCD as for QED (with 
extra 5"^ color factors appearing in all relevant equations); at (7 = 0, QCD is the same free- 
field theory as QED (except for additional color indices); QED in any of the conventional 
relativistic gauges can be obtained by treating the \\ (S^A^)^ as an 'interaction' (as above). 
Taking advantage of these observations, one can set up QCD in the form used above for 
QED. 

As a final preliminary step, we write 



with f;j^ = df,Al - d^Al and C'[A] = -i (2f;^ + gf^^A^Al) [gf^^AlAl], and for subse- 
quent usage, after an integration-by-parts, we note the exact relation 

- ^ /f^ = 1 i-d') Al + \j {d.Alf + j C'[Al (10) 

(In the next few paragraphs, for simphcity, we suppress the fermion/quark variables, which 
will be re- inserted at the end of this discussion). 

In order to select a particular relativistic gauge, one can multiply the 2nd RHS term of 
( ITOl) by A, and include this term as part of the interaction, thereby obtaining the familiar 
QCD GF in the relativistic gauge specified by 

3SD[7]=A/'e^^^'[^^]-e^'^^'^'''^-eH^-°^°^ (11) 

or, after re- arrangement, 

3SDb]=A^e^^-'[^^]-e^^-°^'-^ (12) 

with the determinantal phase factor of (jlj) included in the normalization A/", and a 5"* 
associated with each free-gluon propagator Dc- 

After re- inserting the quark variables, and after re- arrangement, expansion of f|T2l) in 
powers of g clearly generates the conventional Feynman graphs of perturbation theory in 
the gauge C,. It is clear that all choices of A are possible except A = 1, for that choice leads 
to C oo and an ill-defined gluon propagator. This is an unfortunate situation, because 
the choice A = 1 is precisely the one that corresponds to MGI in QCD, as is clear from ( fTOjl . 

But there is a very simple way of re-writing f lT2|) . by replacing the j C'[A\ of that equation 
by the relation given by (fTOj) . 



^ j C'[A] = --J^' + \jAl i-d') A;-'-J {d,A;)\ (13) 

which (continuing to suppress the quark variables momentarily) yields 

3g),j^[j] = A/- e-3/^'-5(i-^)/(^MA;i)%f/A;i(-a^)A« ^ ^ .ef/rD(°>.,_ (^4) 



It is now obvious that the choice A = 1 can be made. It will become clear below that, using 
(|T^ . the functional operations effectively treat gluons as if they were quanta of a 'ghost' 
field. The gluons of the theory, never measurable by themselves, disappear effectively from 
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the exact calculation of every QCD correlation function, without approximation and without 
exception. This 'ghost mechanism' occurs because all factors of 
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of f|T4|) are, in the sum of all virtual gluon processes, effectively removed by the action of 
the term 

of (IT^ . In the end, the gluon acts as a 'spark plug' to generate the MGI and MLC interac- 
tions of the theory, which then take on remarkably simple forms. 

However, if one is interested in the radiative corrections to the gauge-dependent, free- 
field gluon propagator, which corrections are now guaranteed to be gauge- independent, the 
leading RHS factor of 

62 ■ 



should be retained in the rearranged expression of (fT4l) . taken at A = 1, 

Otherwise, and for the specific examples to follow, this factor plays no role and will be 
suppressed. The resulting GF is then MGI, and its superscript will be suppressed. 
After re-inserting quark variables, ( ITSll becomes 



3QCD[j,ry,r/]=Ar6-^/A-D^"'-A.e-t/^^+i/^-H^)-^.6^/^-«^[^]-^+^[^l , (16) 

and we next invoke the representation suggested by Halpern 

e-y^' =M' j d[x]6i^(^-)'+^^^^'^^-, (17) 

where 

i a 

SO that f|T7|) represents a functional integral over the anti-symmetric tensor field x^u- Fol" 
lowing the standard definition [11], all space-time is broken up into small cells of size 5^ 
about each point Wj, and A^' is a normalization constant so chosen such that the right hand 
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side of (fT7j) becomes equal to unity as F^^^ — )■ 0. In this way, the GF may be re- written as 

3QCD[j,^-,r]] =Ar / d[x]ei/-^ .,H!^-^H!H-n^.e^!n-^^y^l^^n^\ ,(19) 

where exp [25^^] is the hnkage operator with S^^'-' = — f / jx' Dc^"* ■ 

Quarks and anti-quarks are treated as stable entities during any scattering, production 
or binding process, which means that relevant functional derivatives with respect to the 
sources r/, f], will bring down factors of Gc{x,y\A), one such factor for each quark or anti- 
quark under discussion. For example, by standard mass-shell amputation, one can then 
pass to the construction of a scattering amplitude of a pair of quarks, or of any number 
of quarks, or of a quark-anti-quark pair. If this scattering is to occur at high-energies 
and small momentum transfer, a convenient and relatively simple eikonal approximation is 
available, derived in detail in Appendix B of Ref. j2o|- This will be a most convenient tool 
in Section IVIIt where we employ the well-known connection between an eikonal function 
dependent upon impact parameter, and an effective potential which is the cause of the 
scattering or production or binding which leads to that eikonal function. 

But it is worth emphasizing that any such, simplifying eikonal approximation is not to be 
confused with the exact functional representations corresponding to scattering, production 
and binding. Because the linkage operator in (fT9!) represents an effective Gaussian functional 
operation upon the A-dependence contained within Gc[v4] and 'L[A\, and because there 
exist Fradkin representations sl of these functionals which are Gaussian in A, the specific 
functional operations required, resulting from well-defined functional derivatives with respect 
to the sources, 77, f/, according to the physical process under consideration, can be performed 
exactly. This produces the sum of all Feynman graphs corresponding to the exchange of 
an infinite number of gluons between quarks and/or anti-quarks, exhibited in terms of the 
Fradkin functional parameters that define the representations for Gc[A] and L[A]. For 
clarity and convenience, we reproduce exact expressions of these Fradkin representations in 
Appendix El 

The result of this approach is gauge-invariant and contains new structures and require- 
ments of a non-perturbative nature. It displays results of such simplicity that one is able, 
as in Section IVIIt to present an analytic derivation of quark binding potentials. A separate 
paper using this approach will present an analytic derivation of a nucleon-nucleon scattering 



and (deuteron) binding potential, which is, to our knowledge, the very first analytic example 
of Nuclear Physics derived from basic QCD. 



II. EXPLICIT GAUGE INVARIANCE 

The correlation functions of QCD are obtained by appropriate functional differentiation 
of f[T^ with respect to gluon and quark sources. Since we are here concerned only with 
quark (Q) or anti-quark (Q) interactions, in which all possible numbers of virtual gluons are 
exchanged, we immediately set the gluon sources j equal to zero. All Q/Q amplitudes are 
then obtained by pairwise functional differentiation of the quark sources ri{y), fi{x), and each 
such operation "brings down" one of a set of (properly anti-symmetrized) Green's functions, 
Gc[v4]. For example, the 2-point quark propagator will involve the FI J d[x] and the linkage 
operator acting upon Gc{x, y\A) exp {L[yl]}, followed by setting A ^ 0. Similarly, the Q/Q 
scattering amplitude will be obtained from the same functional operations acting upon the 
(anti-symmetrized) combination Gc[A]Gc[A] exp {L[A]}, followed by A — )• 0, as 



and other fermionic 2n-point functions are obtained in the same way. 

The Fradkin functional representations for Gc[A] and L[yl], derived in Appendix El dis- 
play a Gaussian dependence on A, and hence the linkage operations of ( l20l) . in any order 
of the expansion of exp{L[y4]| in powers of L, or somewhat more conveniently, using a 
functional cluster expansion j?! for expjS)^''} operating upon exp{L[A]}, can be performed 
exactly. One small complication, easily surmounted, is due to the non-Abelian nature of 
QCD, in which the A-dependence appears inside an ordered exponential, ordered in terms 
of an invariant "Schwinger proper time" variable, s, as in the explicit representation for 
Gc{x, y\A) and ^i[A]. However cumbersome these representations may appear, the essential 
fact is that they are Gaussian in A, and hence the linkage operations, corresponding to the 
summation over all possible gluon exchanges between Q and/or Q lines, may be performed 



M{xi,yi;x2,y2) 



(20) 




A=0 
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exactly. One then finds a simple, effectively local representation for the sum over all such ex- 
changes, here called a "gluon bundle". Henceforth, "bundle diagrams" will replace Feynman 
diagrams containing individually-specified gluon exchanges. 

To display the gauge invariance of all such QCD correlation functions, one first notes 
that L[A] is explicitly invariant under arbitrary changes of the full QCD gauge group Q]. 
Then, it is convenient to combine the Gaussian A-dependence of every entering Gc[^] into 
the quantity 



exp 



d'zQ^,iz)AVz) 



(21) 



where K and Q are local functions of the Fradkin variables, collectively denoted by u^{s'), 
and the i^"(si), ^"(^2), $J^^(si) and $^^(52) are needed to extract the — u{s')) from 
ordered exponentials. Note that the Q and K are also to represent the sum of similar 
contributions from each of the Gc{x,y\A) which collectively generate the amplitude under 
consideration. For example, in the case of the 4-point function one will obtain from the 
product of Gc,i[A] and Gc,i[v4], 

(22) 

"32 



Kitiz) = 2/ / ds'6^'\z -y, + n(s'))r'''=<l'^,,i(.') 

+2/ r ds'S^'\z -y, + tx(.0)r'^$^.,n(^') 
Jo 



and 



q,%{z) = -2gd.^t^^,iz) - g H ds' 5^'\z - y^ + n(s')) u'^is')n;^,is') 

Jo 

-2gd,^:^j{z) - g rds'6^^\z-y, + u{s'))u'(s')nij{s'), 



(23) 



where subscripts 1,2 and I, E are used interchangeably to denote particles 1 and 2; and, for 
clarity, u{s') is used to denote particle 2, and the function 



r ds' S^'\z - y + u{s')) ^';,is') 
Jo 



(24) 



is introduced in Q for ease of presentation. For higher quark n-point functions, there will 
be additional terms contributing to Q and K, but their forms will be the same. 

Combining the quadratic and linear A-dependence from K and Q above, and including 
that L[74] dependence explicitly written in f[T^ . the needed linkage operation reads 



exp 



6 



6 



5A 



5A 



exp 



- / A-K-A + i 



QA 



■exp(L[v4]), 



(25) 
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where 



K?.(^) + gr'^X^iz)] - z') + {z\ (Df)-^ \z') 



(26) 



and 



^l{z) = Q,l{z) + d,xl,{z). 



(27) 



In K, all terms but the inverse of the gluon propagator are local. 

Eq. f l25|) requires the linkage operator to act upon the product of two functional of A, 
and this can be represented by the identity 



(28) 



where, with an obvious notation, the 'cross- linkage' operator exp{ 2) } is defined by 



(29) 



With the identifications 



J^i[A] = exp 



A-K-A + i QA 



, T2[A]=exp{L[A]), 



(30) 



the evaluation of e®-* J-"! [A] is given by a standard, functional identity p, l7| , and reads 



e^^J^i[A] = exp 



^ /" Q-DW-(l-K-DW)-i-Q-iTrln(l-D,-K) 



(31) 



■ exp 



^ /A.K.(l-Df.K)-^A 



q-(i-k.dW)-^a 



It is to be noted that the kernel D^c^ ■ ( 1 — K • D 



»(0) 



reduces to 



dW . (1 - K ■ B^y' = ■ [l - (K + (dW)-^) . D 



-1 



K \ 



(32) 



where instead of fl26p . one now has 



(33) 



It is interesting to note that the additive character of ( 13311 in both spin and 'isospin'-related 
tensors appears as a natural structure of this formulation jisl ]. 
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In the limit A 0, Eq. ([28]) now yields 



(34) 



exp 



- 1 Q ■ K-^ ■ Q + ixr InK + Itt In (-Df) 



2 / M' ' M' 



■ exp 



exp 



Now observe that the first exponential term on the second line of is exactly exp {—Da'}, 
and serves to remove the exp {2}^'} of the operation (exp {Dyi'} • J-2[yl']). With the exception 



of an irrelevant exp 



Tr In ( -D^°) 



factor, to be absorbed into an overall normalization, 



what remains to all orders of coupling for every such process is therefore the generic structure 



e^-^ Ti[A]T2[A] = AT exp 



-- / Q K-^ Q + -Tr InK 



(35) 



■ exp 



2 6A 5A 1^ 5 A 



exp(L[A]), 



in which the now-useless prime of A' has been suppressed. From fl35|) one may now draw 
the following conclusions: 

(i) Nothing in Eq. ( I35l) refers to Dc°'' , which means that gauge invariance is here rigorously 
achieved as a matter of gauge independence. Such invariance cannot be more manifest. The 



importance and nove^ 
for example, of Ref. 



ty of this exact formulation can be recognized and appreciated in view. 



nj and (l2| 



(ii) As expressed in Eq. ( |35|) . each linkage operation upon exp {L[A]} now consists of the 
exchange of a full 'bundle' of gluons, represented by a factor of K^^, while all of the cubic 
and quartic gluon interactions are conveniently incorporated into the remaining Halpern 
functional integral. This is quite different, in structure and interpretation, from linkages 
involving the corresponding L[y4] of an abelian theory, which entail the exchange of but a 
single virtual boson. It now becomes convenient to replace Feynman diagrams, containing 
different orders of gluon exchange, by 'bundle diagrams', in which every graph depicts the 
non-perturbative exchange of all possible gluons between a pair of quark and/or anti-quark 
lines. This will be described, in detail, in the next Section. 

It should be noted that quite similar forms involving at least the part {gf ■ x)~^ of K^^, 
were previously obtained in an instanton approximation to a gauge-dependent functional 
integral over gluon fluctuations The present result, Eq. (|35|1 . shows that such forms are 
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an integral part of the exact QCD theory. This manifest gauge-invariant construction does 
not work for QED, where the simple rearrangement leading from Eq. f|T2|) to f lT4|) cannot be 
implemented. 

(iii) A striking aspect of Eq. fl5^ is that, because K = K.+(gf-x) and the K and Q coming 
from li[A] are all local functions, with non-zero matrix elements (2;|K|2') = K.(z) S^^\z — z'), 
the contributions of Eq. (l35l) will depend on the Fradkin and Halpern variables in a specific 
but local way. This remarkable property will hereafter be called 'Effective Locality' (EL), and 
one can now display the practical usefulness of this description in which relevant functional 
integrals can effectively be reduced to a few sets of ordinary integrals js]. 



III. EFFECTIVE LOCALITY 



It is worth emphasizing the locality aspects of the results in the previous Section. Perhaps, 
the simplest example of EL is in the context of the simple, but non-trivial, quenched eikonal 
scattering model studied in Ref. j4|, which contained the Halpern functional integral 



A/" / d[x] ■ [det {gf ■ x)] ^ ■ exp 



Q-igf-xr'-Q 



(36) 



where the neglect of terms contributing to possible quark self-energies, in the limit of strong 
coupling, replace the exponential factor of (15^ by the argument 



03 ^^I- 



(37) 



The Fradkin -u'-variables have been replaced by asymptotic 4-momenta by virtue of the 
simplifying eikonal approximation, and those momenta are automatically cancelled as will 
be made clear after (1551) and In (157|1 . the color factors ^^"(O), fii(O) are peaked at 

Si = and S2 = 0, as will be noted shortly in conjunction with the function f{b), which 
depends on the collision's impact parameter b. 

In Ref. 4[, as a result of EL, the argument of the Halpern variable xi"^) of (137|) was 
shown to be fixed at a specific value wq = {0,y±,0) [26], where y denotes the CM space- 
time coordinate of one of the scattering quark or antiquark. This is the only xi"^) that is 
relevant to the interaction. All of the other x(^)? for w ^ Wq (and surrounded by a small 
volume of amount 6^, as in the definition/construction of a functional integral ll|]) are 
simply removed from the problem along with their normalization factors, leaving a single, 
normalized functional integral over d"'x('^o)- 
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Further, were the values of y to be subsequently changed, so that — )■ wi, then all the 
w ^ Wi, would become irrelevant, and with their normalization factors will cancel 
away. In effect, the space-time index w is deprived of any physical meaning, and in this sim- 
ple, quenched eikonal scattering, can be omitted [When quenching is removed and different 
factors of L[A] are introduced, there will be more than one such J d"x(w) to be performed; 
but the essential, simplifying effect of EL remains]. In brief, thanks to the EL property of 
this formulation, the Halpern functional integral can be reduced to a set of ordinary finite 
dimensional integrals ^, which can be evaluated numerically, or approximated by a relevant 
physical approximation. 



IV. QCD TRANSVERSE FLUCTUATIONS 

A basic distinction between QCD and other theories, in particular QED, must now be 
made, following from the materials of Sec. [Ill and IIIII It is central to all applications to 
follow. 

That distinction occurs because the quanta of all (abelian) quantized fields may be ex- 
pected to obey standard quantum-mechanical properties, such as perfect position depen- 
dence at the cost of unknown momenta, and vice-versa. But this is impossible for quarks 
since they always appear asymptotically in bound states. Their transverse positions or mo- 
menta can never, in principle, be exactly measured. Neglect of this distinction produces an 
absurdity in the exact evaluation of all QCD amplitudes. 

A phenomenological change in the basic QCD Lagrangian will accordingly be proposed, 
such that a probability amplitude of transverse fluctuations is automatically contained in 
the new Lagrangian, which eventually leads to those potentials essential to quark binding 
into hadrons, and hadron binding into nuclei. Then, all absurdities in estimates of QCD 
amplitudes are removed, enabling one to analytically calculate the effective potentials that 
produce quark binding, as in Sec. IVIIl as well as nucleon scattering and binding potentials 

Before proceeding with this phenomenological change, it is useful to describe the obstruc- 
tion which will occur. As a typical and important part of a 4-point fermionic function, one 
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finds an exponential factor of 

+ '-g j d'^w ^ dsi ^ ds2 <(si)<(s2) (38) 



X - 2/1 + 6^^\w - 1/2 + ^(S2)). 

This expression, corresponding to the interaction of particles 1 and 2, is obtained in the 
approximation of quenching and by neglecting quark's spins: We emphasize that the full 
non-approximate expression dsiplays exactly the same forms as the one under consideration, 
which is why the point can be made using this simplified example. 

In fl38l) . u^, fij, and Si are variables associated with G\{xi,yi\A), whereas u^, and 
S2 refer to G\{x2,y2\A). The last line of (l38ll may be written as 

5^'\w - 2/1 + u{si)) 5^^\yi - 2/2 + ^(^2) - u{si)) (39) 

and one sees that, as a consequence of EL, this interaction is peaked at wq = yi — u{si). 
This means that all the other w ^ wq are irrelevant to the interaction, and, as discussed 
above, are removed along with their normalization factors, leaving dependence only upon a 
functional integration at one single point wq. However, use of the Image Measure Theorem 8| 
can convert this functional integral into an ordinary, finite-dimensional integral, hereinafter 
denoted by / d^xiy^o)- 

The point central to the argument of this Section has now being reached, as one evaluates 
the support of the second delta- function in ( l39l) . which may be expressed as the product of 
delta-functions, in time, longitudinal and transverse coordinates, 

^(l/lO-l/20 + Mo(s2) -Mo(si)) (40) 
X ^{yi-L - 2/2L + Uh{s2) - Ml(si)) 

X 5'^'^\yiL - mi. + uj_{s2) - uj_{si)). 

In the CM of quark 1 and quark 2, one can choose the origin of each time coordinate as 
the time of their closest approach, and then the time difference yio — y2o is always zero. If 
the Qi and Q2 are scattering, then ?/il + ?/2L = 0, since their longitudinal projections are in 
opposite directions; alternatively, if the Qi and Q2 are bound together, then ?/il = ?/2L, and 
their difference vanishes. Either choice makes no difference at all to the following analysis, 
and so we adopt the simplest, second possibihty, ?/il — 2/2L = 0. 
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But then, how should one interpret such a factor as 



5(So(s2) - no(si)) ? (41) 

At face value, at any given couple of values (si, S2) G ]0, s] x]0, s], and any pair of arbitrary 
functions {u, u), each belonging to some infinite dimensional functional space, the probability 
of coincidence of u{si) with u{s2) is likely to be infinitesimally small if not zero. 

Consider the time-coordinate delta-function, 6{uo{s2) — 'Uo(si)), which can have a zero 
argument whenever Uo{s2) and Uo{si) coincide. Assume this happens at a set of points si, 
so that 

5(mo(s2) -Mo(si)) (42) 

= ^6{uo{si) - Mo(si) + (S2 - Si) ■ UoiSi) H ) 

e 

= E^77-^^(^2-3,) 

i \"'Oy^e)\ uo{st)=uo{si) 

In a similar way, the longitudinal delta-function may be evaluated as 



(43) 



ML(Sm)=«L(s2)-i'«L(sf) 



and their product as 



under the restrictions Uo{sm) = uo{si) and ui^{sm) = ui^{si). Now, uq and Uq, and ul and 
Ml are continuous but otherwise completely arbitrary functions: The probability that the 
intersections of 'Uo(si) with ^0(^2) and of 'Ul(si) with ■ul(s2) occur at exactly the same points 
is therefore arbitrarily small. The only place where all four of these continuous functions have 
the same value is at si = S2 = 0, where, by definition of these functions, u^{0) = u^{0) = 0, 
and hence, this pair of delta- functions collapse to the simple product, 

2 V % Si Mo ^2 Mo Si % S2 y 



Note that in Wiener functional space 23], a proof of ( l45l) can be given by the cogent form of 
Theorem jisl. [igI. Note also that because of the factors appearing in the denominator, when 
the large quark's momenta Pi^ and P21, are substituted for the derivatives of the Fradkin's 
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fields u'^{si) and ul{s2), then the automatic suppression of the momenta dependence in (136|) 
becomes obvious since, up to a sign, they compensate for the very same factors in the first 
line of dM]). 

By definition one has 

u^{s)= I ds'u'^is'), M^(0) = 0, (46) 
Jo 

then it is the remaining, transverse delta-function of (HO!) which is now most relevant, the 
term 

5('Hm-m) = 5('H&), (47) 

where h denotes the impact parameter, or transverse distance between the two scattering 
particles. This 5^^^ (6) appears in the exponential of (!35|) . and the question arises as to what 
meaning it can be assigned. Depending on its argument, a delta-function is either zero or 
infinite. In the first case this means that there is no interaction, while the second case means 



that at 6 = 0, one has an infinite phase factor, suggestive of hard disc scattering |25 |. 

The relevant question is therefore why such a delta-function 5^'^^{h) appears at all. The 
answer is that the assumption has earlier been made, in the conventional abelian way, 
that the quark and/or anti-quark may be treated as ordinary particles, whereas asymptotic 
quarks exist only in bound states: Their transverse coordinates cannot, in principle, be 
specified, and there is therefore no reason to retain the conventional (Abelian) practice in 
which such measurement is assumed possible. This is the interpretation that will be adopted 
here, taking the (5'-^^ (6) outcome as a serious warning that some form of quark transverse 
fluctuation is necessary. 

Why does this happens in QCD? Because QCD possesses EL, which conventional abelian 
theories do not. The latter display sums over interconnected propagators, which provide 
a certain vagueness of position, whereas in the exact non-perturbative QCD, as described 
above, one finds the sharp determination of delta-functions corresponding to the EL property, 
and transverse imprecision must therefore be introduced separately, as a fundamental input 
to the theory. 

In Ref . jj| , it was suggested that this difficulty be treated in an ad hoc phenomenological 
way, by replacing (5'-^^ (6) by the smoothly varying, effective Gaussian 

^(6) = (27r)-2 / d^fce*'-'-^. 
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where /x is a mass parameter on the order of the Q-Q bound state (which we shall call 
a "model pion"), although we were able to obtain the conclusions of that paper without 
specifying the precise form of ip{b). In this article, we face this question directly, by first 
developing a formalism in which transverse quark coordinates cannot be specified, and then 
showing how this formalism removes all such absurdities, such as that of the exponential fac- 
tor of 5^^-' (6) above. But it must be emphasized that our prescription is phenomenological, for 
there remains to be shown how such an approach could be derived from a more fundamental, 
operator- field version of QCD, in which transverse fluctuations would occur automatically, 
perhaps in relation to a possible non-commutative geometrical phase of non-perturbative 
QCD. 

We would like to point out the scale change where the integral J d'^wx^i'w) in the expo- 
nent of the Halpern representation f|T7|) is broken up into small cells of volume (5p/i)^, 



Xi, Xi = Xiwi 



Upon rescaling, Xi i^phYXi: re-expressing all interactions in terms of x', there appears 
in (|38|) the factor {6ph)'^^{h). In Ref. \a\ where transverse imprecision was treated in an ad 
hoc way, the size of 5ph was taken to be M~^, where M corresponded to a very large 
energy associated with the eikonal limit. Here, we ask the more physical question of just 
how small that 5 may be chosen in the light of an actual measurement, and we let Quantum 
Physics provide the answer: That contribution to the 5ph corresponding to a time separation 
should be chosen as that corresponding to a (CM) longitudinal coordinate should be 
— while that corresponding to each of the transverse coordinates should be 

and hence the physical volume (5ph)^ is proportional to l/(/i-E)^. An alternate way 
of expressing this is that, starting from arbitrarily small separations in each coordinate, we 
average each Xi variable over a physically meaningful distance, and call that average the Xi 
contained in the volume {5phY. When a scale change {5phYxi = Xi is subsequently made in 
(156|) . a factor of {SphY will appear multiplying ip{b). 



V. A PHENOMENOLOGICAL EXPRESSION OF TRANSVERSE IMPRECISION 

Perhaps the simplest way of introducing transverse fluctuation is to average that part of 
the QCD Lagrangian dealing with the quark-gluon interaction, so that the transverse posi- 
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tion of the color-charge current operator il)'-)^T°"ip{x) should be averaged over a small range 
by means of an initially unspecified distribution. One can also demand the same imprecision 
for the vector current ip 7^ il){x) and scalar density ipiplx), but these extra requirements seem 
to complicate the presentation, to no real advantage, and will not be considered here. 

We emphasize that we have here chosen perhaps the simplest phenomenological way of 
introducing a measure of such transverse imprecision/fluctuation; and that there are quite 
possibly other, more profound methods of obtaining the same objective. In fact, we think 
it most probable that a lack of anti-commutation of quark operator fields at equal times, 
but at differing transverse coordinates could be the basic point requiring attention. In 
the interests of simplicity and clarity, we ask the readers indulgence for postponing that 
particular investigation. 

Instead of the conventional quark-gluon contribution to the Lagrangian density, 

£qg = -V^ [m + 7^ {d, - tgA';^T^)] ^, (48) 

in which all field operators occur at the same space-time point, and for which gauge invari- 
ance under the standard QCD gauge transformations is obvious, we now adopt a local - in 
time and longitudinal position - but non-local in its transverse coordinates replacement, 

£'qg = tgj d^fXf^ -f^)^(x')7M^;:(x)r"^(x'), (49) 

where the transverse imprecision function (TIF) a{x± — x'j_) is a real, symmetric function 
of its arguments, of significant value only for distances on the order of the inverse of the 
pion mass, x'^ = {xo,x'j_,xi^), and v4^(x) is left untouched. In this formulation, rigorous 
local gauge-invariance is suppressed for the underlying quark fields, whose quanta have 
unmeasurable transverse positions, but the hadrons all constructed from these quanta will 
nevertheless be proper singlets under SU(3). 

In fact, one may argue that required gauge invariance is maintained when quark proper- 
ties are averaged over distances below which no measurement of their transverse properties 
can be physically performed. It would be most attractive, if and when a better formal- 
ism is invented, if strict gauge invariance could be maintained for all values of transverse 
coordinates; but that is not a physical requirement, rather a mathematical nicety. 

One notes that in the contribution of (H9|) to its part of the Action operator, J d'^x Cqg, 
the x± and x'j_ coordinates can be interchanged, which yields an equivalent form in which 
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every ^^(a;) of the original (HHIl is replaced by J d'^x'j_ a{x± — x'j_) This interchange 

allows a very simple extraction of all such transverse imprecision, since both delta-functions 
of (15^ will now be replaced by 

j dVi± a{yi± - y[j_) j dV2± ^{mi. - y2±) (50) 

One small complication of this procedure is that such 'primed' transverse coordinates will 
now appear in the arguments of e.g., iu± — > and the fixed position coordinate w± 
must now itself be varied. But the difference \y± — y[_ \ is effectively bounded by where fi 
appears in the definition of (f{b) below, and it turns out that a negligible error is made when 
w± is replaced by y±. More details are given in Appendix [Bl where it is shown that this 
approximation is justified for the subsequent calculations of quark and nucleon bindings. 

The first so modified delta-function of fl50|) defines the argument wi of ^ind we 

again observe that the final output of the Halpern FI will be an ordinary integral j d^x, 
independent of the choice of Wi. The second delta-function of fISUl) now involves the a- 
dependence, generating in place of the 6^'^^{yi± — y2±), the combination 

= J dVla(yi±-fl)a(f2±-yl)- 
Inserting 2-dimensional Fourier transforms of each 

a{y± - yl) = / e^'-^y--y'-^ ~a{k^), (52) 
the combination (15T|) becomes 

/ (C^ ^^^^ e'>^<y^^-y^'^\ (53) 

From its definition, a is real, and hence fl53|) becomes 

l-^,e^'-'\~a{k)f^v{b), (54) 

which provides the definition of (p{b). Note that while no restriction has been placed on 
the form of a other than that it is real and symmetric, turns out to be independent of 
the direction of b, that is, ip{b) = ^(b). In this way, the improper S^'^\b) is automatically 
replaced by (p{b). 
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FIG. 1: A gluon bundle exchanged between quarks I and E 
VI. BUNDLE DIAGRAMS 

In the above example of quark and/or antiquark scattering, where the infinite num- 
ber of exchanged gluons appears to originate and end at a single space-time point on a 
quark/antiquark line, modulo transverse imprecision, it may be helpful to introduce the con- 
cept of an exchanged 'gluon bundle', as in Figured! Because of the four-dimensional delta 
function S^^\y'i — y'2 — u{si) + ■u(s2)), arising from the product of the pair of delta func- 
tions of (ED]), and of the subsequent analysis which produces flHT]) . the transverse separation 
h = yi± — y2± satisfies the probability distribution (IMj) . The argument w in (/ ■ xl"^))^^ is 
given hy w = wi = y[ — u{si) = y'l ^ yi, in virtue of fH5|) . and, as explained in Sec. Illlt the 
Halpern functional integral reduces to a set of ordinary integrals [8] that are represented by 
the Bundle Diagram of Figure [H 

Here, therefore, it is understood that time and longitudinal coordinates of the end-points 
of the bundle are the same, whereas their transverse coordinates, measured vertically in the 
figure, are separated. Bundle diagrams are not Feynman diagrams, but offer perhaps a more 
efficient way of representing the sum over all of the Feynman graphs corresponding to such 
multiple gluon exchange. 

A slightly more complicated expression describes gluon bundles exchanged between any 
two of three quarks, as in Figure [2], where, because of EL, the tu-coordinates of each of the 
if 'X)~^ entering into the appropriate Halpern functional integral are the same, even though 
the transverse coordinates of the three quarks can be quite different. 

In contrast, were a closed quark loop, corresponding to a simple relaxation of the quenched 
approximation, to appear between a pair of quarks, joined to each external quark line by the 
exchange of a gluon bundle, as in Figure |3l there will now be two distinct sets of ordinary 
Halpern integrals to be evaluated. 
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FIG. 2: Gluon bundles exchanged among three quarks. 




y2 

FIG. 3: Gluon bundles joining two quarks via a closed quark loop. 

As will be seen elsewhere, the effective diagram of Figure |3] will provide us with the 
essential features of the Nucleon-Nucleon potential for separation lengths beyond 2 fm |l4|. 
Interestingly also, a comparison of bundle diagrams such as those depicted in Figures [H [2] 
and [3] will be shown to provide an interesting qualitative understanding of the difference 
between QCD and its pure Yang-Mills truncation. 



VII. QUARK BINDING POTENTIAL 

This Section should be viewed as a continuation of that of Ref. jj], in which a new 
formalism for analytic QCD calculations was suggested, and which can serve to identify 
an effective potential as the result of the exchange of multiple virtual gluons between Q's 
and/or Q's. The word 'effective' here represents the type of potential one would expect 
to find were the situation that of a scattering of conventional Abelian particles, moving 
at large relative velocities in their CM, as appropriate in an eikonal context. The word 
'Abelian' is relevant in this QCD situation, for the analysis of Ref. |^ suggests that for large 
impact parameters the scattering of such non- Abelian Q's and/or Q's is effectively coherent; 
only as the impact parameter decreases do progressively stronger color fluctuations appear 
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and decrease the effective coupling, which then becomes the statement of non-perturbative 
asymptotic freedom in the variable conjugate to the impact parameter. For such scattering, 
the eikonal model is, by far, the simplest method of extracting the relevant Physics. 

This and the next Section represent two, simple, phenomenological extensions of Sec. IVl 



However, they are supported by a more rigorous analysis [15|| where, using quenched and 
eikonal approximations, the exact closed forms of such QCD fermionic amplitudes are ob- 
tained. 

We shall be interested in binding, or "restoring" potentials, effective potentials which 
prevent bound particles from remaining at relatively large separations. Since large impact 
parameters are relevant to large, three-dimensional distances, what is needed is, firstly, a 
statement of the relevant formulae of Ref. [J] at large impact parameter; and, secondly, the 
corresponding eikonal function found in that limit. It is then a relatively simple matter to 
determine an effective potential which would produce the same eikonal. 

A two-body eikonal scattering amplitude will have the form 6| 

r(s, t) = ^ j e'"^-' [l - e'""^"''^] , (55) 

where s and t denote the standard Mandelstam variables, s = —(^1+^2)^5 t = — (pi — K)^ — 
(P_ in the CM of QQ, and where X(s, 6) is the eikonal function. Following the arguments of 
J], we give here the simplest derivation of results whose qualitative form describes that of 
a more accurate analysis [l^ . One reaches 

and the normalization is such that for g(p{b) — > 0, exp [iX] — 1. Here, n = — 1, that is 
n = 8 at Nc = 3. All the integrals over 

j (Tai j d"an j d"fii J d"fin (57) 

are properly normalized, and they connect the Qf, dependence of (!56|) with the Af, Aj 
Gell-Mann matrices placed between initial and final state vectors, which define the type of 
amplitude desired, Q-Q scattering or Q-Q scattering, etc. 

As in [4], the 8-dimensional integral over j d^xoa is arranged into the form of a radial 
integration over the magnitude of XoS) now called R, and normalized angular integrations. 
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We suppress the latter, along with all normalized / d"^! / d^flj dependence, and rewrite 
fl56|) in the form 



e 



^^(^'^) = A/-' / dRR'e''^^^^, (58) 







where the normalization of the i?- integral is given by Af''^ = dRR!^ e*"^, and where (g) 
is a shorthand denoting the effective coupling after all angular and color integrations and 
matrix elements have been calculated. In this calculation, only the linear dependence of gip 
is needed, while the non-vanishing matrix elements between singlet color states will define 
Q-Q scattering. The integer power of 3 in fl58l) . of the term = i?^"^, arises from the 
^Tr In factor of (!35|) . contributing R~^ to ( !56|) . As noted above, we wish to compare this 
eikonal function at fairly large impact parameter with that corresponding to the scattering 
of two, interacting particles, and then infer what form of effective potential corresponds to 



that eikonal. We therefore assume the usual formula js] relating an eikonal to a specified 
potential, 

/oo 
dz 1^(6 + zPl) =7(5) X(6), (59) 
-00 

where 7(5) depends upon the CM energy of the scattering particles (and the nature of their 
interaction), and pl is a unit vector in the direction of longitudinal motion. Equation ( l59l) 
is true for non-relativistic and relativistic scattering, in potential theory and in field theory. 
For large b, our eikonal turns out to be weakly dependent on s, and the 6-independent 
contributions do not contribute to V{r). 

To obtain V{r) from a given X(6) is then simply a matter of reversing the usual cal- 
culation: One computes the Fourier transform X(/c_l) of X(6), and then extends k± to 
\k\ = [k'j_ + kj]^ . The three-dimensional Fourier transform of ^{k) is then proportional 
to V{r). Following Sec.V, we assume that X(6), which depends on the form of (p{b), is a 
function of 6^; and although we begin by asking for the form of ip{b) for large b, it can 
be shown that the small-6 contributions to the Fourier transforms are unimportant for the 
large-r behavior of V{r). 

One further point requires consideration, for the above remarks are valid when the process 
involves only the scattering of two particles. But when initial energies are high enough such 
that inelastic particle production occurs, then the potential corresponding to such production 
must have a negative imaginary component (so that the corresponding S-matrix, exp [—il-Lt] , 
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will have a decaying time dependence) , thereby conserving probability as a diminution of the 
final, two-particle state when production becomes possible. This is also true when the two 
initial particles have the possibility of binding, and forming a state which was not initially 
present; the probability of the final scattering state must diminish. In other words, the 
eikonal function one expects to find may always be characterized by a complex potential 
of form Vr, — iVi, where in the present case, Vi is that potential which can bind the Q 
and Q into a pion: In a scattering calculation, that potential must therefore appear with 
a multiplicative factor of —i (symbolically, — > ^(Vr — iVi), and for this calculation of 
binding, iX ^ Vi = Vb). 

Returning to f l55p . it is assumed that g 7^ 0, so that the "1" term of the integrand 
contributes a 5^'^\q) that can be discarded. This amounts to a physical assumption: Q and 
Q are not bound rigidly, with an unchanging impact parameter; rather, there is a permanent, 
if small, "in-and-out" transverse motion, which represents the bound state, and hence q^O. 
Then, 

r(s,t) ~ [ d^b e"^-' e''^^''\ (60) 



where exp [«X(6)] is given by f l58|) . In j4| it was assumed that (p{b) could be a smooth function 
which vanishes for large b, such as ~ exp [— (/i6)^], but it turns out that an appropriate choice 
is 

¥.(6) = y.(0) e-('^^)'"^ (61) 

where ^ is real, positive and small. For any such choice of ^, ip{b) becomes small as fib ^ 1. 

For large b it is then sensible to expand the exp [i{g)ip{b)/ R] term of (1581) . retaining only 
the linear {g)ip dependence. This gives, in place of (158|) . 

POO 2 

Jo 

= l + iKM'{g)ip{b) + 



(62) 



where 

poo ^ 

t,M' = Af' dR R^ e'^ = 20F(-z) i , (63) 
Jo 

and A/"' = — |. Remembering that both sides of f l62|) are to be integrated over J d^fee*'^ '', for 
g 7^ we can proceed to the identification 

e^xW = tK^'{g}if{b), (64) 
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or 



iX{b) = In [ip{b)] + ■ ■ ■ , (65) 

where In [ip] is large and is (effectively) small, and where the dots will be commented 
shortly. 

Choosing ip{b) = ip{0) e~^^''^'^~^\ one obtains 



zX(6) = -(/i6)2+« + -- 



(66) 



It is then convenient to use the integrals 



and 



poo 

/ dxx^'^ sin(a;) = r(/i) sin(^), |Re(/i)| < 1, 
Jo 2 



(67) 



(68) 



in which, except for obvious poles, the Gamma functions are analytic in /i, and can be 
continued to the needed values. With 



(69) 



iX(A;x) = - j d^ifibf+^e"'^ '' 

POO 

= -(27r) / dbJo{k±b)fi^+^b^-^^, 
Jo 



and the use of the doubling formula for Gamma functions, working everything through, one 
finds 

VB{r) = -— fi'^' r^+« ^,^—4t r(-2 - sin(^), (70) 



r(-i-|) 

and for small enough ^, the confining potential 

yB(r)^e/^(H''"^ 



(71) 



which can be compared to the results of several machine groups 17|]. It is remarkable that the 
choice (1611) can be viewed as part of a Levy-flight probability distribution 2l|] that generalizes 
Gaussian ones in consistency with the famous Central Limit Theorem of statistical physics. 
It is interesting to note that at a physical level, such a distribution addresses the issue of 
transverse quark motions, describing them in a way that can seem quite relevant to a confined 
context. Likewise, fl6T]) could suggest a most fascinating connection to a non- commutative 
geometrical aspect of non-perturbative QCD, and this should be examined elsewhere. 
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There remains the question of what values should be assigned to ^, and to the unknown 
parameter /i/mq. These issues will be dealt with in the next Section. Incidentally, one may 
observe that in (1^ . additive constants to In [y:>{b)] bring no contribution to Vb(t'). 

To calculate the corresponding effective restoring potential when one of the three quarks 
contributing to the QQQ bound state is separated from the other two, or, more simply, when 
all three quarks are forced to separate from each other, one may refer to Equation (22) of , 
giving the eikonal amplitude corresponding to Coulomb three-particle scattering. Here, the 
specifically Coulomb parts of this amplitude may be replaced by the single Halpern integral 
which connects the three quarks to each other, as it connected the Q and Q to each other in 
the calculation above. If particle 2 of this formula enters the scattering with zero transverse 
momentum in the rest frame of particles 1 and 3, the amplitude simplifies to 



V 



oik 



12 



$(&12,&32,&13) + 



(72) 



where 613 = 612 — ^32, irrelevant kinematic factors multiplying these integrals have been 
suppressed, and \I'3 denotes the combination 



1 + 



1 - $(612, 00, 00) 



+ 



1 - $(00,632,00) 



+ 



1 — $(00, 00, b 



13 J 



so that 7^''^'^ is a completely "connected" amplitude. 

Writing $ as exp [iX(6i2, 632)], its defining integral may be rewritten as 



POO 

J\f' I dRR^ exp 
Jo 



■ ^2 

+ ^ [(^)l2<^(^12) + (^)32<^(^32) + {9)i3^{bl3)] 



(73) 



(74) 



and one may ask for the form f[74l) takes when one or more of the bij becomes large. 

When any one of the bij becomes large, <f{bij) becomes small, and its exponential term 
may be expanded. The simplest situation is when they all become large, so that an expansion 
of (17^ is relevant. The first non-zero and leading term in such an expansion which contains 
the (f dependence, and hence the 6jj-dependence, of all three quarks contributing to the 
singlet QQQ state is that term for which the expansion yields 



1 + kM' [{9),Mbi2) + {9)32V{b32) + {9)iMbi3)] + ■ • • , (75) 
so that, suppressing all normahzed angular and color integrations, for gl 7^ 0, gjj 7^ 0, 

^ [{g),Mh2) + {9)sMh2) + {9)^Mbi3)] + ■■■■ (76) 
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If particles 1 and 3 are now additionally separated, keeping the distances 612 and 632 essen- 
tially fixed, then 

^MhM ^ ^_^'(^)^^^(5^3) + ... ^ (77) 
iX(6i2,632) ^ In [^(bis)] H , 

and we are effectively back in the "pion" situation, where the large- impact parameter eikonal 
of Q-Q scattering was calculated. Clearly, (1771) suggests that when any two of the three 
quarks which contribute to the singlet "nucleon" are well separated, there results a restoring 
potential V{rij) of the same form as that of (1711) . There are, of course, corrections to this 
confining potential, those which are obvious from the approximations used above, as well as 
those corresponding to spin and angular momentum, which have been completely neglected. 



VIII. ESTIMATION OF THE "MODEL PION" MASS 



This Section contains an estimate of the effects of the above binding potential for the 
simplest case of the "pion", a ground state of the Q-Q system calculated using the simplest 
minimization technique 18|, as well as a second minimization with respect to a ratio of 
terms introduced in this analysis. 

In this simplest, non-relativistic estimation, the Hamiltonian of this system is given by 

n = 2m+ —p'^ + V(r), (78) 
m 

where 1/m denotes the inverse of the reduced mass of this equal-mass quark system. The 
ground-state energy is then given by the replacement of p by 1/r, and the subsequent 
minimization of the eigenvalue of this Hamiltonian with respect to r. Assuming ,^ <^ 1, one 
obtains 

2 



/--(j^l (79) 



and 

3^ [2 



£. = 2™ + -{^^-^j . (80) 

Here, /i and ^ are the two parameters required by this analysis in order to have a non-zero 
binding potential; and while their existence has been introduced as a necessary assumption 
for the interacting Q-Q system, there is another unknown quantity, the quark mass m, which 
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must play an important role. If we take the position that /i and ^ will appear as a result 
of a more fundamental QFT (in which the transverse momenta or position coordinates of 
the quanta of the QCD fields can never be measured with precision), then we are free to 
consider which values of m, or of the ratio fi/m, can lead to the lowest value of E. In so 
doing, one is asking if there might exist a dynamical reason for the relatively small value of 
the pion mass, in this approximation to the actual pion; that is, whether approximate chiral 
symmetry has a dynamical basis. 

Following this approach and using the small-^ limit of (IHOl) . 



mJ V 2 / \m 



(81) 



Then, upon variation with respect to x = fi/m, the function Eo{x) will have an extremum 
/ \ 1/2 

at xq = 2^/'* f I j , and since 



d^Ea(x) 



dx"^ 



2-^/^1) '''(2-^ )>0, (82) 



that extremum is a minimum. Substituting the value oi fx/m = 1?!^ ( j ) into (1811) yields 



Eo = /ie'/'2-i/^[l + 3], (83) 

from which one infers that there is three times as much energy in the gluon field and Q- 
Q-kinetic energies as in the quark rest masses. Intuitively, one expects that E ~ ~ /i, 
which suggests that ^ ~ \/2/\Q. And finally, one then has an estimate of the quark mass, 
in terms of ^ and /i. 

Of course, this double-minimization calculation of Eq cannot be exactly identified with 
the precise pion mass, because of the approximations made above, as well as the omission of 
more complicated singlet terms, such as the contributions coming from QQQQ terms. But 
the Physics seems to be reasonably correct. Note also that we are dealing with "textbook" 
QCD, containing but one type of quark and its complement of massless, SU(3) gluons; flavors 
and electro-weak corrections can be added, as desired. 

A similar minimization analysis can be made for the nucleon ground state, but there then 
are other degrees of freedom which should be included, and treated properly by a serious 
attempt at a solution of such a three-body, relativistic problem. This we must leave to 
others whose numerical abilities far exceed our own. 
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IX. SUMMARY AND SPECULATION 



The previous Sections have presented the formulation of a new approach to analytic, non- 
perturbative, gauge-invariant QCD, from the first step of demonstrating how gauge invari- 
ance of the Schwinger/Symanzik generating functional can be assured; to the explicit demon- 
stration of how that invariance is achieved by gauge-independence of the non-perturbative 
sums over all gluons exchanged between scattering quarks and/or antiquarks; to the expres- 
sion of a new property of non-perturbative QCD called Effective Locality, and the necessity 
of introducing 'transverse imprecision' in the basic Lagrangian of QCD; and ending with 
an explicit determination of the quark-binding potential, and its use in estimating a model 
pion. 

All of this has been accomplished in the context of a 'textbook' QCD example, involving 
but one type of massive quark and its complement of massless SU(3) gluons. For ease 
of presentation, a most convenient eikonal model has been used for high-energy processes, 
although all steps could have been carried out without approximation by the use of an exact 
Fradkin represenation for Gc[A]. Flavors and electro- weak interactions can be included when 
desired. The method of assuring gauge invariance, although simple, had been overlooked for 
decades; while the EL property, along with the necessity of transverse fluctuations, are, to 
the best of our knowledge, new requirements inherent to non-perturbative QCD. 

As an indication of the usefulness of this approach, a subsequent paper now under prepa- 
ration will derive from 'realistic' QCD, a quite good representation of a nucleon-nucleon 
binding potential, as appropriate to the deuteron. This, to our knowledge at least, is the 
first analytic QCD derivation relevant to the entire subject of Nuclear Physics, while it 
also suggests obvious glueball candidates. That calculation requires the use of a closed- 
quark- loop, as in Fig. |3l and the subject of further 'radiative corrections', and indeed, of 
non-perturbative QCD renormalization, makes clear that the adjective 'non-perturbative' 
used in this paper refers only to the complete summation of all possible gluon exchanges 
between quarks and/or anti-quarks. The questions then arises: Is it possible to extend this 
formalism to include sums over all closed quark loops? What is the effect of gluon bundle 
exchange within and between closed quark loops? Is it possible to achieve a totally non- 
perturbative statement of any given QCD amplitude? We hope to answer these and other, 
related questions in subsequent publications. 
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Appendix A: Fradkin's Representations for Gc[^] and L[A] 



The exact functional representations of these two functionals of A{x) are perhaps the most 
useful tools in all of QFT, for they allow that A-dependence of these functionals to be ex- 
tracted from inside ordered exponentials; and because they, themselves, are Gaussian in their 
dependence upon A{x), they permit the functional operations of the Schwinger/Symanzik 
generating functional (Gaussian functional integration, or functional linkage operation) to 
be performed exactly. This corresponds to an explicit sum over all Feynman graphs rele- 
vant to the process under consideration, with the results expressed in terms of functional 
integrals over the Fradkin variables; and in the present QCD case, because of EL, those 
non-perturbative results can be extracted and related to physical measurements. 

The causal quark Green's function (which is essentially the most customary Feynman 
one) can be written as js, 6| 



Gc[A] = [m + i7 ■ n][m + (7 ■ Uf]-^ = [m + i-f- U] ■ i 



ds e"*"'"" e^"(^-") 



where 11 = i[dfj_ — iQ-^r''] and (7 ■ U)^ = + iga^y F^^^r" with a^y 
Fradkin's method p, |6| and replacing with i-^, one obtains 



(Al) 

WliJi^lA- Following 



Gc{x,y\A) 

ds e-^"""' 



( (,1ds'—i 



m - 7/, 



S{x-y+ / ds' v{s')) 



X exp 



-ig I ds' 
'0 



v^is')Aliy- I v)T^^ia^yYlXy- 



(A2) 









1). 


/o 
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Then, one can insert a functional 'resolution of unity' of form 

1= / d[u]6[u{s')- r ds" v{s")], 



(A3) 



and replace the delta-functional 6[u{s') — ds" v{s")] with a functional integral over 



and then the Green's function becomes 



19| 



Ge(x,i/|A) (A4) 



where h{si,S2) = ds' Q{si — s')Q{s2 — s'). To remove the A-dependence out of the lin- 
ear (mass) term, one can replace igA'^(y — u{s))t'^ with —5^27^ operating on the ordered 
exponential so that 



Gc{x,y\A) 

POO 

i ds e~ 
Jo 



(A5) 



X 



6 



To extract the A-dependence out of the ordered exponential, one may use the following 
identities, 



1 = I d[a]S[a'^is')+gu'^is')A;{y-uis'))]., 
1 = / d[S]6[S;,{s')-gF;,{y-u{s'))], 



(A6) 



and the ordered exponential becomes 



-iff ds' u'^is') A-iy-uis')) T^+g ds'a^, F", ^ 



(A7) 



^ / ds' is') a- {s')~i J ds' is') 3-,(s') 

^ ^-ig J ds' u'^(s') n-is') A-iy-uis'))+ig J ds' F-Jy-u(s')) 



where N'n and A/$ are constants that normalize the functional representations of the delta- 
functionals. All A-dependence is removed from the ordered exponential and the resulting 
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form of the Green's function is exact (it entails no approximation). Alternatively, extracting 
the A-dependence out of the ordered exponential can also be achieved by using the functional 
translation operator, and one writes 



H-5.0 



For the closed- fermion-loop functional L[A], one can write [6| 
L[A]- 



1 






2 







{^ = 0}}, 



(A8) 



(A9) 



where the trace TV sums over all degrees of freedom, space-time coordinates, spin and color. 
The Fradkin representation proceeds along the same steps as in the case of Gc[v4], and the 
closed-fermion-loop functional reads 

1 r°° ds 

s 



L[A] 



g -ism^ g-|Trln{2/i) 



(AlO) 



-{9 = 0}, 

where the trace tr sums over color and spinor indices. Also, Fradkin's variables have been 
denoted by v{s'), instead of u{s'), in order to distinguish them from those appearing in the 
Green's function Gc[A]. One finds 

1 /" ds ^_ism'2 ^-^Trln(2fe) 



L[A] = — / _e— 6-2^^1^^/^) (All) 
2 Jo -5 

xMnM^ j d^x J d[a] j dp] j d[H] j d[$] 

X j rfH5W(t;(.))et/o'^^'[^'(^')]' 

^ ^-i ! ds' n-(s')a-{s')-i J ds' ^l,{s')=.l,{s') .^^ f^^i ds' [a-{s')-ia^,,S-^(s')\r-^ 
Xe-*S !a ds' v'^is') n-is') A-(x-v(s'))-2igJ d'^z (9.*S^(^)) A-(z) 



+ig^ J ds' rb-^l,(s') Ab^(x-v(s')) Al{x-v{s')) 



xe 
{9 = 0} 



where the same properties as those of Gc[A] can be read off readily. 
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Appendix B: Effective Locality versus Transverse Average 



Before transverse imprecision was introduced, EL had the effect of attaching to the rep- 
resentative symbol [/ ■ of each gluon bundle, exchanged between quark and/or 
antiquark of respective CM coordinates yi and y2, a pair of delta functions, S^^\w — yi + 
u{si))5^'^\yi — y2 + u{s2) — u{si)) , as used in the text, or the pair 5^^\w — yi + sipi)6^^\yi — 
2/2 + S2P2 — S1P2) as used in an eikonal approximation {4]. For either case one finds fixed 
values of wq and wl, and w± = yi± = —y2±- Then, as claimed in the text, the Halpern 
FI can be reduced to an ordinary set of J d"x integrals. In the process, though, one makes 
a systematic error, of the eikonal-type, by neglecting variations of the impact parameter 
or, correspondingly, of momentum transfer in the core parts of the matrix element. In the 
context of the exact expression of the first pair of delta functions above, that ad hoc approx- 
imation avoided the much more complicated analysis of the transverse Fradkin's difference 

With transverse imprecision now being included, the situation changes for the better, 
in the sense that no such approximation need be made. But this change now requires a 
slightly more complicated justification of the argument which replaces Halpern's FI by a set 
of ordinary integrals. For the question arises if this useful simplification is also true when 
the ■w± inside the [/ ■ xi"^)]^^ factor is itself given by y^, and is being integrated over the 
J d'^y' in that exponential factor, as in the discussion of the text leading to (153|) . It was 
there noted that the replacement of that w± by yi± or —y2± is a reasonable approximation. 

The following argument is intended to give that approximation a more detailed justifica- 
tion. 

Consider the Halpern FI 

AT / d[x] [det (/ ■ xT^' 6^^'^""'^'+^^^'^'^" "(^^i^-^")"(^^2x-27l)[/-x«)]-^^ (Bl) 



where y'^ = {yo]yL,y±)y the normalization is defined so that the FI of f IBip equals 1 when 
g = 0. The dependence of color, time and longitudinal coordinate has been omitted for 
simplification of presentation. 

As in the definition of this or any such FI, f d^wx"^ is understood to mean 6^ J2f=i 
where the subscript i denotes the value of x at the space-time point w±i, and 6^ corresponds 
to a small volume surrounding that point, which is to become arbitrarily small as becomes 



35 



arbitrarily large As mentioned in the text by the end of Sec. |Vl residual (5-dependence 
will be re-expressed in terms of physically significant quantities as a last step; but for the 
following argument, all the transverse coordinate differences are to be taken as arbitrarily 
small. 

Now, re-scale the variables such that = Xe, and re-write (IBip as 



Let us also break up the J d'^y'i integral into an infinite series of terms: one is free to choose 
the individual y'^ coordinates as exactly those which define the transverse positions of the 
Xe = xi'^i)- III this way, ( ]B2|) may be re- written as 

J\f / d[x] [det (/ ■ x)]~^ T.exj+igS'^ T,t<yi±-y'Le)^iy-2±-y±e)[f-xiy'Le)]''^^ 



where A^^ is understood as a true infinitesimal quantity, and where, for simplicity, we 
suppress explicit dependence on yo and But now (]B3P may be written as the product of 
integrals, 

^ - r 

N 

^ JI ¥{ig5^ AJ^ a{y,^ - y'^,) a{y,^ - y[,)) 
I 

N 

e 

where flB4p denotes the normalized product of all such (y^^)-valued integrals, and ¥{zi) 
denotes the ordinary integral J d^Xi over the variable associated with y\_^. That integral is 
well defined in the sense that, for l^;^! < 1, as is the case here, it can be expressed as an 
absolutely-convergent series, or as a converging integral over a set of eigenvalues in a random 
matrix calculation. 

One then expects to be able to write F(z) in terms of its Fourier transform, 

¥{ze) = J dg¥{g)e'''', (B5) 

where the normalization condition of (]B4p stipulates that / dg¥{g) = 1. Since Zi is propor- 
tional to the infinitesimal A^^, one may expand in powers of Z£, 



¥{ze) = I dg¥{g) [1 + izeg + ■■■] = 1 + i / d^ ^F(^) ze + ■ ■ ■ , (B6) 
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so that (lB4p becomes approximately 



n 



1 + i dg g¥{g) ze 



(B7) 



l + i I Aggf{g) ^ 



l + i ^F(^) -igS d y'^ a{yi^ - y[) a{y2± - v'l) 



With 



dVl - y±) ^iy2± - y±) 

Eq. (]B4I) becomes 



|a(g)| e 



2 Jqiyi±-y2±) 



Vib), (B8) 



(B9) 



which is just the first-order expansion of the result obtained in the text when y'^ was shifted 
to ylj_ or —y2±- And since gS"^ ip is expected to be small, 6'^ ip 1, it is, in effect, equivalent 
to 

j dg¥{g) e^d^aS^^^^b) ^ 

which is just the integral of (]B2p when the intuitively equivalent change y'^ — )■ ylj_ has been 
made in the argument of (/ ■ and after the residual 6^ dependence has been continued 
to the measurably-significant value of l/{fiE). Only the first-order form, corresponding to 
( ]B9|) is needed in the calculations of quark- and nucleon-binding. 



[1] H. Reinhardt, K. Langfeld, and L. v. Smekal, Phys. Lett. B300, 111 (1993). 

[2] R. Hofmann, Int. J. Mod. Phys. A20, 4123 (2006); Erratum-ibid. A21 6515 (2006). 

[3] D. D. Ferrante, G. S. Guralnik, Z. Guralnik and C. Pehlevan, Complex Path Integrals and the 

Space of Theories, Brown Univsersity Preprmt: BROWN-HET-1611 (2011). 
[4] H. M. Fried, Y. Gabellini, T. Grandou and Y.-M. Sheu, Eur. Phys. J. C65, 395 (2010). 
[5] E. S. Fradkin, Nucl. Phys. 76, 588 (1966). 

[6] H. M. Fried, Basics of Functional Methods and Eikonal Models (Editions Prontieres, Gif-sur- 

Yvette Cedex, France 1990) 
[7] H. M. Fried, Functional Methods and Models in Quantum Field Theory (The MIT Press, 

Cambridge, MA 1972) 



37 



[8] B. Candelpergher and T. Grandou, work in progress. 

[9] M. B. Halpern, Phys. Rev. D16, 1798 (1977). 

[10] M. B. Halpern, Phys. Rev. D16, 3515 (1977). 

[11] A. Zee, Quantum Field Theory in a Nutshell (Princeton University Press, Princeton 2010) 

[12] A. Guay, Geometrical apects of local gauge symmetry, (2004), 

jhttp : //philsci-archive . pitt . edu/id/eprint72T33[ 

[13] K. Huang and D. R. Stump, Phys. Rev. Lett. 37, 545 (1976); Phys. Rev. D15, 3660 (1977). 



[14] H. M. Fried, Y. Gabelhrn, T. Grandou and Y.-M. Sheu, arXiv:1203.6137v l [hep-ph] (2012). 
[15] H. M. Fried, T. Grandou and Y.-M. Sheu, work in completion. 

[16] T. Grandou, On Some Aspects of the QCD Effective Locality, Proceedings of the Eleventh 

Workshop on Non-Perturbative QCD, Paris, June 2011, Edited by B. Miiller and C. I. Tan, 

|http : //www . slac . Stanford . edu/econf /C1106064/1 
[17] See, for example, the seminal papers by Y. Nambu, Phys. Lett. B80, 372 (1979); M. Luscher, 

K. Symanzik, P. Weisz, Nucl. Phys. B173, 365 (1980); M. Luscher, Nucl. Phys. B180, 317 

(1981). 

[18] F. Balibar, A. Laverne and J. M. Levy Leblond, Quantique: Elements, 
|http : //eel . archives-ouvertes . f r/docs/00/13/61/89/PDF/elem_5f ev07 . pdf [ 

[19] Y.-M. Sheu, Finite- Temperature Quantum Electrodynamics: General Theory and Bloch- 
Nordsieck Estimates of Fermion Damping in a Hot Medium, PhD Thesis, Brown University, 
May 2008. 

[20] H. M. Fried, Y. Gabellini and J. Avan, Eur. Phys. J. C13, 699 (2000). 

[21] For example, see http : //en . wikipedia . org/ wiki/Stable_distributioii} 



[22] H. M. Fried, Green's Functions and Ordered Exponentials (Cambridge University Press, Cam- 
bridge 2002) 

[23] G. W. Johnson and M. L. Lapidus, The Feynman Integral and Feynman's Operational Calculus 

(Oxford University Press, Oxford 2000) 
[24] L S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (Academic Press, 

London 1994) formula 8.253.1 
[25] L. I. Schiff, Quantum Mechanics, 3rd Ed. (McGraw-Hih, 1968) 
[26] wq = (Oo,yj_,OL), correcting the expression given in Ref. 
[27] H.M. Fried, K. Kang and B. H. J. McKellar, Phys. Rev. A28, 738 (1983). 



38 



